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Abstract 

Simple, linear equations relate microscopic swimmers to the cor- 
responding gliders and pumps. They have the following set of con- 
sequences: The swimming velocity of free swimmers can be inferred 
from the force on the tethered swimmer and vice versa; A tethered 
swimmer dissipates more energy than a free swimmer; It is possible 
to swim with arbitrarily high efficiency, but it is impossible to pump 
with arbitrarily high efficiency and finally that pumping is geometric. 
We also solve several optimization problems associated with swimming 
and pumping: The problem of optimal anchoring for a certain class 
of swimmers that includes the Purcell swimmer and the three linked 
spheres and the optimal geometries of helices considered as swimmers 
and pumps. 

1 Introduction 

Low Reynolds numbers hydrodynamics governs the locomotion of tiny nat- 
ural swimmers such as micro-organisms and tiny artificial swimmers, such 
as microbots. It also governs micro-pumps [1] , tethered swimmers [2], [3] and 
gliding under the action of external forces. Our purpose here is to discuss 
some interesting consequences of simple, yet general and exact relations that 
relate pumping, swimming and gliding at low Reynolds numbers. The equa- 
tions are not new, and appear for example in [U [5], where they serve to 
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analyze swimming at low Reynolds numbers. The interpretation of these 
equations as a relation between three apparently different objects: Swim- 
mers, pumps and gliders, gives them a new and different flavor and suggests 
new applications which we explore below. 

Our first tool is Eq. ([1]) which relates the velocity of a free swimmer to 
the (torque and) force needed to anchor it to a fixed location. This equation 
appears in the analysis of swimmers [HE]. It can be applied to experiments 
on tethered swimmers [6l [71 [8] where it can be used to infer the velocity of 
a free swimmers from measurements of the force on tethered one (and vice 
versa) and when both are measured, to an estimate of the gliding resistance 
matrix M of Eq. 

Our second tool is Eq. ([7|) which says that the power needed to operate 
a pump — a tethered swimmer — is the sum of the power invested by the 
corresponding free swimmer, and the power needed to tow it. (The result 
holds for arbitrary container of any shape or size). It is useful in studying 
and comparing the efficiencies of swimmers and pumps. In particular, it says 
that a tethered swimmer spends more power than a free swimmer. As we 
shall see it also implies that it is impossible to pump with arbitrary high 
efficiency. 

Physically, both relations may be viewed as an expression of the elemen- 
tary observation that swimmers and pumps are the flip sides of the same 
object: A tethered swimmer is a pump [2] and an unbolted pump swims. 
This observation plays limited role in high Reynolds numbers, because it 
does not lead to any useful equations for the non-linear Navier Stokes equa- 
tions. For microscopic objects, where the equations reduce to the linear 
Stokes equations, the observation translates to linear relations. 

Formally, a pump is a swimmer with a distinguished point — the point of 
anchoring. The pump may have different properties depending on the point 
of anchoring, (we shall see an example where it pumps in different directions 
depending on the anchoring point). We shall solve a problem, posed by E. 
Yariv [3], of how to find the optimal anchoring point for a certain class of 
swimmers, which we call "linear swimmers". The class includes the Purcel 
swimmer [H], the "three linked spheres" [TU] and pushmepuUyou [TTj . 

There are various notions of optimality that one needs to consider in 
the study of swimmers and pumps. For swimmers, one notion of optimality 
is to maximize the distance covered in one stroke. A different notion is 
to minimize the dissipation in covering a given distance at a given speed. 
Similarly, for a pump one may want to maximize the momentum transfer to 
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the fluid in one stroke, or, alternatively, to minimize the dissipation for given 
total momentum transfer and rate. How are optimal swimmers related to 
optimal pumps? We study this question for helices and show that optimal 
swimmers and pumps have different geometry: There are four optimal pitch 
angles depending on what one optimizes and whether the helix is viewed as 
a swimmer or a pump. 

2 Anchoring as a choice of gauge 

Let us start by briefly discussing the gauge issues that arise when considering 
swimming [12], pumping and gliding. 

Fixing a gauge is already an issue for a glider. A glider is a rigid body 
undergoing an Euclidean motion under the action of an external force. There 
is no canonical way to decompose a general Euclidean motion into a trans- 
lation and a rotation [13] . Such a decomposition requires choosing a flducial 
reference point in the body to flx the translation. In the theory of rigid body 
a natural choice is the center of mass. This is, however, not a natural choice 
at low Reynolds number. This is because this is the limit when inertia plays 
no role and the glider may be viewed as massless. Since the choice becomes 
arbitrary, it may be viewed as a choice of gauge. 

In the case of a swimmer one needs to flx a point arbitrarily and in 
addition, to flx a flducial frame. This is because a swimmer is a deformable 
body and such a frame is required to flx the rotation [12]. This is, again, a 
choice of gauge. 

A pump is an anchored swimmer with a distinguished point and a dis- 
tinguished frame which are determined by the way the pump is anchored. 
When we write equation that involve swimmers, pumps and gliders, we pick 
the flducial point and frame determine by the the way the pump is anchored. 

3 Triahty of swimmer pumps and ghders 

We shall flrst recall the derivation of the linear relation between the 6 di- 
mensional force-torque vector, Fp = [Fp, Np) which keeps a pump anchored 
with flxed position and orientation, and the 6 dimensional velocity-angular- 
velocity vector, Vs = {Vs,uJs) associated with the corresponding autonomous 
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swimmer. 

Fp = -MV„ (1) 

This holds for all times {Fp,M and are time dependent quantities). M 
is a 6 X 6 matrix of linear-transport coefficients of the corresponding glider. 

F, = MV, (2) 

namely, the corresponding rigid body, moving at (generalized) velocity Vg 
under the action of the (generalized) force F^. Note the change in sign. The 
matrix M depends on the geometry of the body. It is a positive matrix of 
the form [Ti] : 

where K, C, and fl are 3x3 real matrices. Note that Eq. ([2]) fails in two 
dimensions |15j . 

Eq. (JT]) follows from the linearity of the Stokes equations and the no-slip 
boundary conditions. Let (9S denote the surface of the device. Any vector 
field on 9S can be decomposed into a deformation and a rigid body motion 
as follows: Any rigid motion is of the form = V + x x. Pick V to be 
the velocity of a fiducial point and u the rotation of the fiducial frame. The 
deformation field is then, by definition, what remains when the rigid motion 
is subtracted from the given field v. 

Now, decompose the velocity field on the surface of a swimmer, to 
a deformation and rigid-motion as above. The deformation field can be 
identified with the velocity field at the surface of the corresponding pump Vp 
since the pump is anchored with the fiducial frame that neither moves nor 
rotates. The remaining rigid motion is then naturally identified with the 
velocity field on the surface of the glider. The three vector fields are then 
related by 

Vs = Vp -F- Vg, Vg = Vs + ujs X ^ (4) 

where Vg is (by definition) the swimming velocity and Ug the velocity of 
rotation. 

Each of the three velocity fields on (9S, (plus the no-slip zero boundary 
conditions on the surface of the container, if there is one), uniquely determine 
the corresponding velocity field and pressure {v,p) throughout the fiuid. The 
stress tensor, Hij, depends linearly on {v,p) [15] . 
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By the linearity of the Stokes, djiTij = 0, and incompressibihty djVj = 
equations, it is clear that Vg = Vg + Vp and ps = Pg + Pp and then also 

= vTp + Tig. Since = /^^ HijdSj is the drag force acting on the device we 
get that the three force vectors are also linearly related: Fg = Fp + Fg, and 
similarly for the torques. This is summarized by the force-torque identity 
Fg = Fp + Fg. Since the force and torque on an autonomous Stokes swimmer 
vanish, Eq. ([T]) follows from Eq. ([2]). 

Eq. has the following consequences: 

• Micro-Pumping and Micro-Stirring is geometric: The momentum and 
angular momentum transfer in a cycle of a pump, J Fpdt, is indepen- 
dent of its (time) parametrization. In particular, it is independent of 
how fast the pump runs. This is because swimming is geometric [H [12] 
and the matrix M is a function of the pumping cycle, but not of its 
parametrization. 

• Scallop theorem for pumps: One can not swim at low Reynolds numbers 
with self- retracing strokes. This is known as the "Scallop theorem" [S] . 
An analog for pumps states that there is neither momentum nor angular 
momentum transfer in a pumping cycle that is self-retracing. This can 
be seen from the fact that dt is balanced by — dt when the path 
is retraced, and this remains true for M dt. 

We shall now derive an equation, originally due to [5], which relates the 
power expenditure of swimmers, pumps and gliders. It follows from Lorentz 
reciprocity for Stokes flows that says that if {vj,7ijk) and (w^-,7r^-^) are the 
velocity and stress fields for two solutions of the Stokes equations in the 
domain S then [T3|: 

/ V'i TTij dSj = Vi TT-j dSj (5) 

JdT. JdT. 

For the problem at hand, we may take 9S to be the surface of our device 
(since the velocity fields vanish on the rest of the boundary associated with 
the container). The area element dS is chosen normal to the surface and 
pointing into the fluid. Now apply the Lorentz reciprocity to a pump and a 
swimmer velocity fields and use Eq. (j4]) on both sides. This gives 

-Pg + Yg-Fg = -Pp -Yg-Fp (6) 
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where P, is the power invested by the swimmer and Pp the power invested 
by the pump. Since the force and torque on the swimmer vanish, Fg = 0, we 
get, using Eq. ([1]) a hnear relation between the powers: 

Pp-Ps = -V, ■ Fp = V, ■ MVs = Pc,>0 (7) 

Pg is the power needed to tow the ghder. Since both swimming and towing 
require positive power, at any moment pumping is more costly than swim- 
ming or dragging. 

The linearity of Eq. ([7]) is a noteworthy, and somewhat unexpected. 
Eq. (jlj) says that the corresponding velocity fields are linearly related. Since 
power at low Reynolds number is quadratic in the velocity a linear relation 
between the powers is not what one may naively expect. 

The most interesting consequence of this relation which, at least for us, 
was somewhat of a surprise, is that a pumps needs more power than a swim- 
mer and so the power consumed by a tethered swimmer is actually an upper 
bound on the power consumed by it when freed. 

One of the remarkable facts about low Reynolds number swimming is 
that even though the dynamics is governed by dissipation it is still possible 
to swim with arbitrarily high efficiency [HI [16]. For this to happen, the 
swimming velocity should be non-zero, but the energy dissipation of the 
swimmer, Pg, should be zero, thus Pp = Pg > 0. The same cannot happened 
for pumps: for a pump to be with arbitrarily hight efficiency, it must have 
non-zero momentum transfer to the fluid and zero power. Since M of Eq. ([3]) 
is a strictly positive matrix Pg is quadratic in the force and can not vanish 
if the pump transfers momentum to the fluid. 

4 Linear swimmers and optimal anchoring 

E. Yariv [1] posed the following problem: Find the anchoring point which 
optimizes the momentum transfer to the fluid. The general case is compli- 
cated. A class of swimmer for which this question can be answered relatively 
easily is the class of "linear swimmers" — swimmers made of segments, so 
that the velocity of different points on the same segment depend linearly on 
the distance between the points. This class contains the three linked spheres 
[TU], the pushmepuUyou [TT], Purcell's three linked swimmer [TTl [TB], the 
'N-Linked' swimmer [T^, but not, for example, the treadmiller [TB]. 
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Figure 1: A treadmiller transports its skin from head to back. The tread- 
miller in the figure swims to the left while the motion of its skin guarantees 
that the ambient fiuid is left almost undisturbed. If it is anchored at the 
gray point then it transfers momentum to the fiuid (to the right). The power 
of towing a frozen treadmiller and of pumping almost coincide refiecting the 
fact that a treadmiller swims with little dissipation. 



Shifting the anchoring point by r, the resistance matrixes K, C and Vt of 
Eq. E] will be changed ioKr = K,Cr = C- KR and VLr = VL- RKR + CR- 
RC^ [llj, where R is the antisymmetric matrix associated with the vector 
r. The change in Vg is linear in r, and Ug is unchanged [13]. Using Eq. ([1]), 
the change in the force Fp (and hence in the linear-momentum transfer to 
the fiuid) is linear with r, and the optimal point which maximized the force 
must be at one of the edges of the segments. 

(The change in the torque Tp (and hence in the angular-momentum trans- 
fer to the fiuid) is quadratic in f, so the maximum can be either at the edges, 
or at the point in which the differential of the angular momentum with re- 
spect to the anchoring point is zero (in cases where there is such a point). In 
either case, one has to check only a few points.) 

A case in point is the "three linked spheres" [TU]. There are three can- 
didates for the optimal anchoring point: The three spheres. The two outer 
spheres are related by symmetry so the two interesting cases are either an- 
choring on an outer sphere or in the middle sphere. Detailed calculations 
show [20] that the maximum momentum transfer will be when the anchoring 
point is on one of the outer spheres. When the swimmer is anchored on 
the middle sphere, the momentum transfer turns out to be in the opposite 
direction. Calculation for the efficiency shows that the optimal point in this 
case is also at the outer spheres. 
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5 Helices as swimmers and pumps 



Eq. ([T]) and Eq. ([7]) have an analogs for non-autonomous rigid swimmers such 
as a hehx rotating by the action of an external torque [21]. This is a case 
that is very easy to treat separately. From Eq. ([2]) applied to the helix twice, 
once as swimmer and once as a pump we get the analog of Eq. ([T]) : 

Fp = Cuj = -KVs (8) 

The analog of Eq. ([7]) follows immediately from the definition of the power 
P = -F ■ V and Eq. ^ again: 

Pp - = Fp ■ Vp - F, ■ V, = -Vs ■Fp = Pg (9) 

It follows from this that the difference in power between a swimmer and a 
pump is minimized, for given swimming velocity, if the swimming direction 
coincides with the smallest eigenvalue of K which is the direction of optimal 
gliding. 

The distinction between optimal pumps and optimal swimmers [22] can 
be nicely illustrated by the considering the example of a rotating helix, or a 
screw. This motion has been studied carefully in [231 [21] whose analysis goes 
beyond what we need here. For a thin helix the slender-body theory of Cox 
[25] disposes of much of the hard work. Cox theory has the small parameter 
(log/t)"^ where k is typically the ratio of the (large) radius of curvature of 
the slender body, r in the case of a helix, to its (small) diameter. To leading 
order in (log the local force field on the body is fixed by the local velocity 
field: 

dF{x) = k{t{t ■ v) - 2^r)dx, A; = ^ (10) 

t{x) is a unit tangent vector to the slender-body at x and v(a;) the velocity 
of the point x of the body. This result may be interpreted as the statement 
that each line element has twice the resistance in the transverse direction 
than the resistance in the longitudinal direction, and that the motion of one 
line element does not affect the force on another element (to leading order). 

Consider a helix of radius r, pitch angle 6 and total length £. The helix 
is described by the parameterized curve 

(r cos(/), r sin0, tsin^), (j) = -cos6, tG[0,i] (11) 
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Suppose the helix is being rotated at frequency uj about its axis. Substituting 
the velocity field of a rotating helix, with an unknown swimming velocity in 
the z-direction, into Eq. ffTOj) . and setting the total force in the z-direction 
to zero, fixes the swimming velocity. Dotting the force with the velocity and 
integrating gives the power. This slightly tedious calculation gives for the 
swimming velocity (along the axis) and the power of swimming: 



K sin 2^ 



ijjr 3 + cos IQ kiu'^r^ 3 + cos 29 
Similarly, for the pumping force and power one finds 



(12) 



^ = s,n«cos«. ^ = l+sin'« (13) 

Eq. fll2p and fll3l) have the following consequences for optimizing pumps and 
swimmers: 

• Given ur, a swimmer velocity Vg is maximized at pitch angle 6 = 54.74°. 

• Given ur, the pumping force Fp is maximized at 6* = 45°. 

Consider now optimizing both the pitch angle 6 and rotation frequency uj so 
that the swimming velocity is maximized for a given power. Namely 

max{K I Ps = const} (14) 

and similarly for pumping, except that Fp replaces Vs and Pp replaces P^. 
A simple calculation shows that this is equivalent to optimizing V^/Ps and 
Fp / Pp with respect to 9. (These ratios are independent of uj and so invariant 
under scaling time). One then finds: 

• The efficiency of swimming, V^/Ps, is optimized at ^ = 49.9°. The 
efficiency is proportional to {k€)~^ which favors small swimmers in less 
viscous media, as one physically expects. 

• The efficiency of pumping, Fp/Pp, is optimized at ^ = 42.9°. The 
efficiency is proportional to {k€) which favors big pumps at more viscous 
media. Micro-pumps are perforce inefficient. 
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Figure 2: A Purcell three linked swimmer, left, controls the two angles 6i and 
62- When bolted it effectively splits into two independent wind-shield wipers 
each of which is self retracing. It will, therefore, not pump. PushmepuUyou, 
right, controls the distance between the two spheres and the ratio of their 
volumes. It can have arbitrarily large swimming efficiency. 



There is a somewhat unrelated, yet insightful fact that one learns from 
the above computation regarding the difference between motion in a very 
viscous fluid and motion in a solid. The naive intuition that the two are 
similar at very high viscosity would imply that a helix moves like a cork- 
screw and so would move one pitch in one turn. This is actually never the 
case, no matter how large fi is. In fact, the ratio of velocities of a helix to a 
cork-screw is independent of fi and by Eq. (IT^ 

= ^Q^^ < 1 (15) 
ur sin 6 1 + cos^ 6 ^ 

A helix needs at least two turns to advance the distance of its threads. 



6 Pumps that do no swim 

We have noted that the "three linked spheres" can pump either to the right 
or to the left depending on whether it is anchored on the center sphere or on 
an external sphere. There is, therefore, an intermediate point where it will 
not pump. 

There are also pumps that will not swim: Evidently, if the swimming 
stroke is right-left and up-down symmetric, the swimmer will not move by 
symmetry. It can, however, be bolted in a way that breaks the symmetry 
to give an effective pump. For example - consider a right-left symmetric 
pushmepuUyou jlT], where the two spheres inflate and deflate in phase. By 
symmetry, it will not swim, however, anchoring in any point except the 
middle point will lead to net momentum transfer. 
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